Magnetic field-induced insulating behavior in highly oriented pyrolitic graphite 



D. V. Khveshchenko 

Department of Physics and Astronomy, University of North Carolina, Chapel Hill, NC 27599 

We propose an explanation for the apparent semimetal-insulator transition observed in highly 
oriented pyrolitic graphite in the presence of magnetic field perpendicular to the layers. We show 
that the magnetic field opens an excitonic gap in the linear spectrum of the Coulomb interacting 
quasiparticles, in a close analogy with the phenomenon of dynamical chiral symmetry breaking in 
the relativistic theories of the 2 + 1-dimensional Dirac fermions. Our strong-coupling appoach allows 
for a non-perturbative description of the corresponding critical behavior. 
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The recently discovered carbon-based materials pro- 
vide a new playground for applications of the advanced 
methods of quantum field theory. Although the best 
known example is that of the one-dimensional carbon 
nanotubes described as the Luttinger liquid, some of the 
available non-perturbative techniques can also be applied 
to higher dimensional systems, such as layered highly ori- 
ented pyrolitic graphite (HOPG). 

In a single sheet of graphite, the low-energy spec- 
trum of the quasiparticle excitations becomes linear in 
the vicinity of the two conical points located at the 
opposite corners of the two-dimensional Brillouin zone 
where the conduction and valence bands touch upon each 
other These low-energy excitations can be conve- 
niently described in terms of a four-component Dirac 
spinor = (ipiAv, ipiBo, ^2Aa, ip2B<r), combining the 
Bloch states ^io-(r) with spin a which are composed of 
the momenta near one of the conical points (i — 1,2) 
and propagate independently on the two different sub- 
lattices (r = A, B) of the bi-partite hexagonal lattice of 
the graphite sheet. In the following discussion, we will 
treat the number of the spin components N as a tunable 
parameter, the physical case corresponding to N = 2. 

The use of the Dirac spinor representation allows one 
to cast the quasiparticle kinetic energy in the relativistic- 
like form 
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where ^> a — ^tlo- The reducible representation of 
the 4x4 7-matrices 70,1,2 = (73, ir-2, — iri) ® T3 sat- 
isfying the usual anticommutation relations {7^,7^} = 
2diag(l, — 1, — 1) is given in terms of the triplet of the 
Pauli matrices r it and the velocity v ~ 10 6 m/s is pro- 
portional to the width of the electronic 7r-orbital band. 

The Lorentz invariance of the non-interacting Hamil- 
tonian is not, however, respected by the interaction term 
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which accounts for the long-ranged part of the Coulomb 
coupling whose strength is characterized by the dimen- 
sionless parameter g = 2ire 2 /cqv. 

The earlier perturbative studies of the effects of the 
Coulomb interaction resulted in the prediction that, upon 
renormalization, the strength of the effective coupling 
g(e) ~ l/|lne| monotonically decreases with the energy 
e, hence the paramagnetic semimetallic ground state re- 
mains stable 0]. However, this conclusion appears to 
contradict the recent experimental observation of a fer- 
romagnetic bulk magnetization (inconsistent with the es- 
timated number of magnetic impurities) in the HOPG 
samples showing insulating behavior of the resistivity [^) . 

In fact, the large value of the bare Coulomb coupling 
g ft 10 suggests that a more appropriate starting point 
might the strong coupling regime where perturbation the- 
ory fails and a more capable approach is needed. In 
the present paper, we propose such an approach by fo- 
cussing on a recent experimental observation of the ap- 
parent magnetic field-driven semimetal-insulator transi- 
tion in HOPG 0j and demonstrate that external mag- 
netic field can trigger the instability towards excitonic 
insulator phase. Interestingly enough, the latter appears 
to have much in common with the phenomenon of chi- 
ral symmetry breaking (CSB) which has been previously 
studied in the relativistic theories of the interacting Dirac 
fermions. 

In the case of graphite, the issue of CSB comes about 
due to the invariance of Eqs.(l,2) under arbitrary U(2N) 
rotations of the 2iV-component vector comprised of the 
chiral Dirac fermions $ (l,r)o- = |(1± 75)^0-, where the 
matrix 75 = 1 ® T2 anticommutes with any 7 M . 

In a quantum system, strong interactions can give rise 
to the appearance of a fermion mass and gapping of the 
fcrmion spectrum, thereby breaking the continuous chi- 
ral symmetry U(2N) down to its subgroup U (N) ®U(N) 
which corresponds to the independent rotations of ^ La 
and ^ Ra . 

As one important 2 + 1-dimensional example, CSB can 
be caused by the interaction with a scalar Higgs- Yukawa 
(HY) bosonic mode coupled to the Dirac fermions via the 
mass operator ^a^a- I n t ms case, CSB is known to 
occur for any number of fermion species N, provided that 
the strength of the (intrinsically attractive) HY coupling 
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exceeds a certain critical value. Recently, this model has 
been applied to the analysis of the condensed matter sys- 
tems where the scalar bosonic field describes fluctuations 
of a superconducting order parameter || or piesoelectric 
phonons 

In contrast, the repulsive Lorentz-invariant vector-like 
coupling via the current operator \E f 'y Ai \l/ drives the 
Dirac fermions towards the CSB transition, regardless 
of the coupling strength, provided that the number of 
fermion species N is sufficiently small (N < N c ). This 
behavior is believed to occur in the strong coupling in- 
frared fixed point in the 2+1-dimcnsional Quantum Elec- 
trodynamics (QED3) where the zero temperature value 
of N c was found to be smaller or equal to 3/2 0. 

However, this situation changes drastically in the pres- 
ence of magnetic field which suppresses the orbital mo- 
tion of the Dirac fermions and collapses their spectrum 
into a discrete set of the (relativistic) Landau levels. This 
effectively reduces the dimensionality of the problem, 
thus enabling CSB to occur, regardless of the coupling 
strength and/or the number N of fermion species. 

In the Dirac picture of the quasiparticle excitations in 
layered graphite, the onset of CSB would be manifested 
by a non-zero value of the order parameter ^f„^„ = 

YtiA^iAa^iA" ~ TiBtr^iBa)- lt determines the magni- 
tude of the fermion gap proportional to the electron den- 
sity imbalance between the A and B sublattices and cor- 
responds to the formation of a site-centered charge den- 
sity wave (CDW) in the excitonic insulating ground state. 

While for a sufficiently small N the excitonic instabil- 
ity could develop even in a single layer of graphite || , in 
the physical case N — 2 it is unlikely to occur in the ab- 
sence of the inter-layer Coulomb repulsion. Indeed, in a 
realistic HOPG system consisting of many layers stacked 
in a staggered configuration, the latter favors sponta- 
neous depletion of the electron density on one of the two 
sublattices (e.g., A) formed by the carbon atoms posi- 
tioned, vertically, at the centers (respectively, corners) of 
the hexagons in the adjacent layers. Within each layer, 
such a depletion (and excess occupation of the comple- 
mentary sublattice, e.g., B) conforms to one of the two 
degenerate CDW patterns which then alternate between 
the layers, thereby keeping the electrons in the adjacent 
layers as far apart as possible and further strengthen- 
ing the propensity towards the excitonic instability. Al- 
though the minimal strength of the inter-layer Coulomb 
repulsion required for CSB to occur in HOPG remains 
unknown, this whole situation, too, changes when the 
system is exposed to magnetic field normal to the layers 
which promotes CSB even in the absence of the inter- 
layer coupling. 

In the presence of the magnetic field B, the Dirac 
fermion Green function remains diagonal in the space 
of the physical spin (not to be confused with the 4 ® 4 
space of the 7-matrices representing the orbital dynam- 



ics), and, for a given a, it reads as 
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where the translationally non-invariant phase factor con- 
tains the vector potential of the external field A fl (x) — 
(0, —Bxijl, Bxi/2). Upon separating this factor out, 
one obtains a translationally (albeit, not Lorentz-) in- 
variant Green function Q{x). 

The effect of the fermion interactions can be fully ac- 
counted for by introducing the gap function A(p) as well 
as the wave function (Z) and velocity (Z v ) renormaliza- 
tion factors into the Fourier transform of Q (x) given by 
the integral representation (hereafter, we use the units 
v — e = h = 1 and the relativistic notations, such as 
= ( e >P)) 
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which the exact Green function naturally inherits from 
the bare one, Go{p), given by Eq.(4) with Z = Z v = 1 
and A(p) = 0. 

Owing to its non-perturbative nature, the phenomenon 
of CSB eludes weak-coupling analysis based on perturba- 
tion theory Nonetheless, akin its relativistic counterpart, 
the occurrence of CSB in the system of the Coulomb 
interacting Dirac fermions can be inferred from the non- 
perturbative solution of the Dyson equation for the renor- 
malized Green function 



g- 1 ( J p)-Gv 1 (p)=ig 



d 3 k z ioG(p + k)% 
(2tt) 3 \k\+Ng X (k) 



(5) 



In Eq.(5) we made use of the Ward identity relating the 
vertex function JqZ to the energy derivative of (? _1 (p) 
and cast the effective intra-layer Coulomb interaction in 
the form governed by the scalar (density-density) com- 
ponent of the fermion polarization operator 
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As a result of the broken Lorentz invariance due to both, 
the non-relativistic nature of the Coulomb interaction 
and the presence of the magnetic field, the solution of 
the gap equation (5) can feature totally different depen- 
dencies on the energy and momentum variables. Apart 
from the Zeeman shift g\xbB of the position of the Fermi 
level for the spin-cr fermions, Eq.(5) remains spin degen- 
erate. 

The analysis of Eq.(5) is complicated by the fact that, 
unlike in the previous studies of the excitonic transi- 
tion in semimetals with overlapping conduction and va- 
lence bands, the naive picture of static Debye screening 
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(Xk ~ const) fails to properly describe the feedback of the 
planar Dirac fermions on the bare Coulomb interaction. 
Instead, in the undoped or lightly doped graphite with 
a low density of carriers the zero temperature fermion 
polarization (6) can only be expressed in the form of a 
cumbersome double integral 
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Nonetheless, a progress towards obtaining the solution of 
Eq.(5) can still be made in the strong field limit where 
the distance betweent the adjacent Landau levels of the 
non-interacting Dirac fermions E n = ±yj2\n\B by far ex- 
ceeds the Coulomb interaction-related energy gap A, and 
the only relevant fermion states appear to be those of the 
so-called lowest Landau level (LLL) with n — 0. In this 
regard, the problem at hand bears a certain resemblance 
to the Fractional Quantum Hall effect (FQHE) in the 
system of non-relativistic fermions with a parabolic dis- 
persion. In contrast to the spatially homogeneous FQHE, 
however, the sought solution of the gap equation corre- 
sponds to a non-uniform CDW ground state, as explained 
above. 

In the strong field approximation A <C \^B the fermion 
Green function reduces to its LLL projection 
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and, correspondingly, the zero temperature fermion po- 
larization receives its main contribution from the transi- 
tions between the LLL and the first excited Landau level 
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Neglecting the wave function, velocity, and vertex renor- 
malizations in the scalar part of Eq.(5) (see below) one 
readily obtains a closed equation for the gap function 
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Eq.(10) should be contrasted with that derived in the 
strong field limit in the case of QED 3 where the LLL 
projection eliminates all but the scalar component of the 
Lorentz-invariant abelian gauge interaction, thus result- 
ing in the gap equation (10) with the term |k| in the 
denominator of the integrand replaced by k 2 = k 2 — uj 2 , 



as required by the Lorentz invariance of the bare gauge 
field spectrum. 

It is this difference which gives rise to the non-trivial 
energy dependence of the solution of the QED$ gap equa- 
tion found in Ref. M. By contrast, in our case the solu- 
tion of Eq.(10) remains independent of the energy vari- 
able as long as e < \AB and, being a sole function of the 
momentum, it falls off faster than e~ 2p "/ s for |p| > \[B. 
It can be also shown that the approximate constancy 
of A(p) at small momenta and energies justifies our ne- 
glecting the wave function and velocity renormalization 
{Z = Z V = 1). 

The magnitude of the zero-momentum gap A = A(0) 
estimated from Eq.(10) 
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demonstrates that the strong field condition holds at 
large ./V regardless of the field strength B (which should, 
however, be smaller than the width of the linear part of 
the bare quasiparticle spectrum). In the physical case 
N = 2, a numerical analysis which does not rely on 
the LLL approximation Q shows a qualitatively simi- 
lar behavior, thus supporting the conclusions drawn in 
the large- N limit. 

At finite temperatures the gap decreases, and its mo- 
mentum dependence at |p| < (T^fB) 1 / 2 becomes even 
less pronounced. Although the prohibitive form of the 
finite temperature fermion polarization impedes analyti- 
cal calculations, its main effect can be taken into account 
through the modified lower cutoff in the logarithmically 
divergent momentum integral in Eq.(10) which now be- 
comes the larger of \^B/gN and (Ty/B) 1 / 2 . 

After being extended to finite temperatures, Eq.(10) 
yields the self-consistent equation for the magnitude of 
the gap at small momenta 



A, 



In- 



B At 

_ — tanh— (12) 

4tt7V maxiVB/gN^TVB) 1 / 2 ) 2T 



Eq.(12) possesses a non-trivial solution below the transi- 
tion temperature T C (B) whose large- N estimate is given 
by the expression 
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Conversely, Eq.(13) determines a threshold magnetic 
field B C (T) cx T 2 which has to be exceeded in order for 
CSB to occur at a non-zero temperature T. This relation 
defines a critical line in the B — T phase diagram, along 
which the gap vanishes as 



A T (B -> B c ) cx y/B - B C (T) 



(14) 



As far as the nature of the CSB transition is concerned, 
the dependence (14) is characteristic of the second order 
transition, whereas the conjectured zero field transition 
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(for N < N c ) is of the topological (Kosterlitz-Thouless) 
kind §. 

These predictions should be compared with the avail- 
able experimental evidence obtained from the HOPG 
samples showing metallic behavior of the zero field re- 
sistivity. The data of Ref. indicate that, albeit ab- 
sent in zero field, the apparent semimetal-insulator tran- 
sition can be induced by magnetic field normal to the 
layers. The insulating behavior was found to set in at 
the field-dependent characteristic temperature fitted as 
T* oc \/B — Bo which, apart from the offset field Bq, 
agrees with Eq.(13) and falls into the experimental range 
< lOOif for the applied magnetic field B < 0.2T (in order 
to avoid confusion we remind that in our quasi-relativistic 
system the role of the "speed of light" is played by v). 
The proposed orbital (as opposed to the spin-related) na- 
ture of the observed phenomenon is also consistent with 
the findings of Ref. Q], according to which the in-plane 
magnetic field has a substantially (two orders of magni- 
tude) weaker effect. 

In fact, the spin degeneracy between the triplet and 
singlet excitonic gaps (for N = 2) will be lifted upon in- 
cluding the short-ranged Coulomb exchange interaction 
omitted in Eq.(2) which involves transitions between the 
conduction and valence bands. Alongside with the Zee- 
man coupling, the latter favors the triplet excitonic or- 
der parameter, in accord with the Hund's rule. In the 
doped system, the spin up- and down- states will then 
be occupied asymmetrically at low temperatures, result- 
ing in the occurrence of a ferromagnetic spin polariza- 
tion M = Tr(aG) in a window of the electron chemical 
potential /x set by the gap A. The induced ferromag- 
netic moment will then be proportional to /i ~ A as 
well, which compares favorably with the electron spin 
resonance (ESR) data in HOPG showing the presence of 
the magnetization M oc yJB — Bo on the insulator side 
of the field-induced transition E(J . 

Before concluding, we comment on the previous at- 
tempts to apply the alternate scenario of the magnetic 
field-driven CSB in the Lorentz-invariant HY model to 
the analysis of the quasiparticle transport in the mixed 
state of the planar d-wave superconductors |0]. Specif- 
ically, the authors of Ref. |ll]] focussed on the experi- 
mentally observed kink-like feature in the magnetic field 
dependence of the total (inclusive of both the normal 
quasiparticle and phonon contributions) thermal conduc- 
tivity, the position of the kink scaling with temperature 
as B*(T) oc T 2 [jl2| (this behavior was not seen, however, 
in a more recent experiment fig] ). 

On the theoretical side, the suggestion of using the 
standard HY model to describe the nodal ci-wave quasi- 
particles in the mixed state must be taken cautiously. 
Apart from the default choice of the quasiparticle inter- 
action in the form of the attractive HY coupling, the 
"magnetic catalysis" scenario of Refs. |ll]] requires that 
the total effective magnetic field, which the nodal quasi- 



particles are exposed to, has a non-zero overall flux, re- 
sulting in the formation of the Landau levels. 

However, as pointed out by several authors 113] , the 
correct picture of the d-wave quasiparticle spectrum is 
rather that of the extended energy bands, owing to the 
fact that in the mixed state the quasiparticles are ac- 
tually experiencing both, the external magnetic field and 
the solenoidal superfluid flow created by the vortices. On 
average, the two fluxes exactly cancel each other [ O , 
thus rendering inapplicable the standard HY mechanism 
of the " magnetic catalysis" even in the regime of weak- 
to-moderate magnetic fields where the individual vortices 
strongly overlap and the physical magnetic field is almost 
uniform. In the light of the above, HOPG might be the 
only presently known example of a condensed matter sys- 
tem where the phenomenon of " magnetic catalysis" can 
indeed occur. 

In summary, we study the problem of the Coulomb 
interaction-driven electronic instabilities in layered 
graphite in the presence of magnetic field. Elaborat- 
ing on the relativistic-like description of the low-energy 
quasiparticle excitations in a single sheet of graphite, we 
propose a possible explanation for the recently discovered 
field-induced semimetal-insulator transition by showing 
that applied magnetic field induces the excitonic insula- 
tor phase, thus gapping up the quasiparticle spectrum 
and creating a site-centered CDW. 

Since the phenomenon in question was only observed 
in perpendicular magnetic field, we believe that the ex- 
perimental findings of Refs. PJl0[ reveal some novel, in- 
trinsically two-dimensional, physics whose origin is differ- 
ent from that of the (semi)metal-insulator transitions ob- 
served in other carbon-based materials, such as the one- 
dimensional ((Ru,Cs)Cqo) as well as three-dimensional 
(KCqo) alkali doped fullerides. The latter are likely to 
be associated with the structural instabilities and/or ac- 
companied by the onset of the antiferromagnetic spin- 
Peierls state. In order to decisively discriminate between 
these possibilities, such experimental techniques as ESR, 
NMR, x-ray diffraction, and electron photoemission, can 
be further implemented. 

The author is grateful to Y. Kopelevich for communi- 
cating his experimental results prior to publication and 
S. Washburn for a valuable discussion. This research was 
supported by the NSF under Grant No. DMR-0071362. 



[1] G. Semenoff, Phys. Rev. Lett. 53, 2449 (1984); F. D. M. 

Haldane, ibid 61, 2015 (1988). 
[2] J. Gonzalez et al, Nucl. Phys. B424, 595 (1994); Phys. 

Rev. Lett. 77, 3589 (1996); Phys. Rev. B59, 2474 (1999). 
[3] Y. Kopelevich et al, J. Low Temp. Phys. 119, 691 (2000). 



4 



[4] Y. Kopelevich et al, Phys. Solid State 41, 1959 (1999); [9] 

H. Kempa et al, Solid State Commun. 115, 539 (2000). [10] 

[5] D. V. Khveshchenko and J. Paaske, Phys. Rev. Lett. 86, [11] 
4672 (2001). 

[6] A. H. Castro-Neto, Phys. Rev. Lett. 86, 4382 (2001). [12] 

[7] T. Appelquist et al, Phys. Rev. D60, 045003 (1999). [13] 

[8] K. Farakos et al, Phys. Rev. D61, 045005 (2000); J. [14] 

Alexandre et al, ibid D62, 105017 (2000); ibid D63, 

065015 (2001). 



D. V. Khveshchenko and H. Leal, unpu blished 



:ond-mat/0106232 



M. S. Sercheli et al, 
G. W. Semenoffet al, Mod. Phys. Lett. A13, 1143 (1998); 
W. V. Liu, Nucl. Phys. B556, 563 (1999). 
K. Krishana et al, Science 277, 83 (1997). 
Y. Ando et al, Phys. Rev. B62, 626 (2000). 
M. Franz and Z. Tesanovic, Phys. Rev. Lett. 84, 554 
(2000); L. Marinelli et al, Phys. Rev. B62, 3448 (2000); 
O. Vafek et al, ibid B63, 134509 (2001). 



5 



